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Abstract 

A new generalization of the t-J model with a nearest-neigbor hop- 
ping is formulated and solved exactly by the Bethe-ansatz method. 
The model describes the dynamics of spin-S fermions with isotropic 
or anisotropic interactions. In the case S=l the magnetic interaction 
is biquadratic in the spin operators. In contrast to the SU (N) gen- 
eralization of the t-J model, studied previously in the literature, the 
present model possesses beyond a massless excitation also a massive 
one. The physical properties indicate the existence of Cooper-type 
pairs with finite binding energy. 
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The t-J model has emerged as a paradigm for studying the low-energy 
electronic properties of the copper-oxide-based high-temperature supercon- 
ductors [1,2]. Although high-T c cuprates are at least two-dimensional systems 
the one-dimensional version of the model and its generalizations are also in- 
tensively studied since in this case exact results can be derived [3-12]. The t-J 
model describes the dynamics of spin-^ fermionic particles with Hamiltonian 
given by 

H = -f £P (ct,c,- +1> „ + 4+1,^) P + J E (44-1 - n 3 n 3+1 /A) , (1) 
j> j 

— * - 

where Cj >cr is the standard fermion creation operator, Sj = ^dj is the particle- 
spin operator and rij is the particle-number operator at site j. The projection 
operator P excludes the double occupation at each site. Unfortunately the 
exact integrability of (1) is obtained only at the supersymmetric point J = 2t 
[3-6]. At this point the model has no gap and the critical exponents governing 
the long-distance behavior of correlation functions were calculated [7]. These 
results show that for any density of holes the spin-spin correlation functions 
dominate the superconducting ones, and as a consequence the model has 
no superconducting properties. These results were extended to the SU(N) 
generalization of the t — J model of fermions of arbitrary spin S [8 - 10]. 
The integrability of an anisotropic generalization of the SU(N) supersym- 
metric t — J model has been shown [13, 14] and the critical exponents of the 
correlation functions have been calculated [15, 16]. 

In this letter we present a new set of models of strong-correlated electrons 
which are exactly solvable. The first example of these models is the spin-1 
biquadratic t-J model with Hamiltonian given by 



H = -tJ2 P { c l* c 3+i,° + c t+i,<T C 3,<r) P ~ J J2 (Sj-Sj+i) 2 ~ n j n j+i 



3,cr 



(2) 



where now Sj = (Sj, Sj , Sj) are spin-1 Pauli operators located at site j. We 
show that this model is exactly integrable at the special point t = J. Actually 
the above Hamiltonian is the isotropic version of a new family of anisotropic 
models describing the dynamics of spin-S* fermions with Hamiltonian 
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L S 

H = -EE p (i c 3 +M+ii,A S ) p 

j=l s=-S 



L 



5 

2 u a utcl s c j>t cf +1 _ s c j+1 _ t - (1 + e 1 )cosh 7 n i n j+1 

s,t=-S 



(3) 



where L is the lattice size, e,e± = ±1 and the parameters u s , which play 
the role of anisotropies should satisfy u s = l/u_ s (s = —S, —S + 1, . . . , S) 
and 2cosh7 = u 2 _ s + u 2 _ s+1 + ... + u|. The particular case S — \ and 
e = —Si = 1 is the anisotropic version of the supersymmetric t-J model. 
The biquadratic t-J model, at t=J, given in (2) is obtained by choosing in 
(3) S — 1, e — —Ei = 1 and u_i — Uq — U\ — 1. For general spin S the 
magnetic interactions can be written as a polynomial of degree 2S in the spin 
operators. 

The exact integrability of these models, from a mathematical point of 
view, comes from the fact that the Hamiltonian density in (3) is related to 
the generators of Hecke algebras [17], with deformation parameter q given by 
the relation q + l/q — 2 cosh 7. 

The eigenstates and eigenvalues of Hamiltonian (3) can be obtained ex- 
actly within the framework of the Bethe-ansatz method [18-21]. The struc- 
ture of the Bethe-ansatz equations follows from the solution of the two- 
particle problem. The two-electron wave function can be written as a prod- 
uct of two factors: a coordinate wave function (referring to the positions and 
momenta of the particles) and a spin part, the global wave function being 
antisymmetric under the exchange of two particles. The scattering matrix 
can be written in the following form 

Sft,(Ai-A 2 ) = [l + (l + ei)cosh 7 $(Ai-A 2 )]<J ai/3 ^ 

- E 1 U a U / 3>^(\ 1 - \ 2 )5 a -f35f3 / - a >, (4) 

where 

sin A 

$ A = -— — — 5 

sm(A — -ry) 

and Aj (j = 1, 2, n) are suitable particle rapidities related to the momenta 
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{kj} of the electrons by 



% = {"na (A d 7) ' ££l = l\ ( 6 ) 
I -@( A i^7), eei = +1, 

with the function defined by 

0(A;7) = 2arctan(cot7 • tan A) ; -7r < 6(A, 7) < it. (7) 

A necessary and sufficient condition for the applicability of the Bethe-ansatz 
method is the Yang-Baxter equation [18,21]. In our case the S-matrix sat- 
isfies these equations in the non-deformed and q-deformed cases [17]. The 
isotropic case corresponds for S > \ to the g-deformed case where u s = l(s = 
—S, . . . , S) and q + l/q = 2S+1. The underlying Hecke algebra of the model 
implies that differently from the supersymmetric t-J model we should have 
gapped spin excitations for S > 1. Up to our knowledge this model is the 
first example of integrable model with the S-matrix of the form (4) which is 
connected with the Hecke algebra. The Hamiltonian (3) is diagonalized by 
standard procedure by imposing periodic boundary conditions on the Bethe 
function. These boundary conditions can be expressed in terms of the trans- 
fer matrix of the non-uniform model which can be constructed on the basis 
of the S-matrix (4) by using the quantum method of the inverse problem [22, 
23]. The rapidities {Xj} that define a n-particle wave function are obtained 
by solving the equations 



sinh(Aj — ry/2) 



1) A(\j), (8) 



sinh(Aj + ry/2) 
where A (A) is the eigenvalue of the transfer matrix 

n 

^7 } } (A) = En^; +i (^-A), (/wi=a). (9) 

{A}*=1 

It is simple to verify that besides the number of particles n, the magnetiza- 
tion J2j Sj an d the number of paired electrons m are conserved quantities in 
the Hamiltonian (3). Two electrons are paired if they are consecutive elec- 
trons with opposite spins and have no unpaired electron between them. The 
complete diagonalization of the transfer matrix (9) is not a simple problem 
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even in the simplest case S — 1, n — L (see, for example, [24]). It is not diffi- 
cult to convince ourselves that in the interesting physical situation where we 
have low density of holes the ground-state will belong to the sector where we 
have zero magnetization and only pairs of electrons. In this sector m = n/2 
and the diagonalization of the transfer matrix of the inhomogeneous model 
(9) gives for e 1 = — 1, the following equations 



sin(A j + i 7 /2) nL 
sin(Aj — ry/2) 



l)m -i fr sin(Aj - A a + ry/2) 
) ^sin^-A^-ry^)' 



" sin(A J - A a - yy/2) _ ™ sin(A Q - + yy) 
/J sin(A i - A Q + i 7 /2) ~ ^ sin(A a - A^ - ry) ' 

In the case e\ = +1 the first set of equations in (10) should be replaced by 



sin(Aj + ry/2) 
sin(A i - i7/2) 



>_ 1)m -i A sin(A, - A; + yy) ™ sin(A j - A Q - i 7 /2) 
sin(A, - A, - ij) £ sin ( Aj - A a + i 7 /2) 



The total energy and momentum of the model are given in terms of the 
particle rapidities \j in the following form 

n n / sinh^ 7 \ 

E — — 2 cos kj = 2eei cosh 7 — , 

pi 3 pA cosh 7- cos 2A J 

n 

P = (11) 

The equations (10) and (11) have the same structure as those appearing in the 
anisotropic t-J model [15,16] provided a suitable definition of the parameter 
7 is given. It means that in spite of the physical processes in the models 
with S — \ and S > \ being quite different there is a "weak equivalence" 
in Baxter's sense [25] between models with different values of spin S in the 
sector where m = n/2. Of course in the general case this equivalence does 
not exist. 

Although the models are exactly integrable for both signs of e and S\ 
in (3) let us now restrict to the more physically interesting case e — 1 and 
si = —1, where we have attraction among pairs. In this case the ground state 
contains m = n/2 bound pairs characterized by a pair of complex electron 
rapidities 
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A± = ^(u a ±i 7 ), v a = 2A Q . (12) 

The second set of equation in (10) is fullfilled within exponential accuracy 
whereas the first set can be treated in the similar way as in [15,16]. Inserting 
(12) in the first set of equations in (10) and introducing the density function 
p(v) for the distribution of v a in the thermodynamic limit, we obtain the 
linear integral equation 

2np(v) = 7) - J $(u - v'; -f)p(v')dv' (13) 



where 



$f • ) - sinn 2 7 (14) 
cosh 27 — cosh f 

In order to minimize the ground-state energy 

—— = —2e / [2 cosh 7 — sinh7$(f ; 7)] p(v)dv (15) 

Lj J I 

the integration interval I in (13) and (15) has to be chosen symmetrically 
around ir(I = [v , 2-k — vq}. The parameter v is determined by the subsidiary 
condition for the total density p = 2m/ L of electrons 

fp(v)dv=^p. (16) 

To study the superconducting properties of the model under consider- 
ation we calculate the long-distance behaviour of the correlation functions 
by finite size studies and application of conformal field theory (see [26-28] 
and references therein). The results of this calculation are the following. 
The long-distance behavior of the density- density and the superconducting 
correlation functions are given by 



(p(r)p(O)) ~ p 2 + A^ 1 + A 2 r- a cos(2k F r); 2k F = irp; (17) 
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p(r) = ]Tc+c r7 , 

7 

G p (r) = (c+ C +. 1) _ 7 c , (5 c 1 ,_ (5 ) ~ Br-?. (18) 



The exponents a and /3 describing the algebraic decay are calculated from 
the dressed charge function £ (v) which is given by the solution of the integral 
equation 

= 1 - ^ / *(v - v'; ^(v')dv', (19) 

and is given by 

a = /T 1 = 2[£M 2 . (20) 

In our one-dimensional system we have no superconductivity in the literal 
sense, since the model does not have finite off-diagonal long-range order. But 
we may say that in our model there is tendency to the superconductivity 
since the superconducting correlations have a longer range than the density- 
density correlations. This happens when (3 < a. Analytically we find a = 2 
for (p = 0) and a — \ for (p = p max — 1)- This implies that for all nonzero 
values of the parameters 7 there is a density regime [0, p c ] where the system 
has dominating superconducting correlations. An analogous behaviour of 
correlation functions can also be observed in the SU(N) generalization of 
the anisotropic t-J model where superconducting properties are caused by 
the introduction of anisotropy in the interactions. However unlike these 
models the superconducting properties in the Hamiltonians (3) are caused 
by both effects, the anisotropy and the value of the spin S (see definition 
of the parameter 7 (3)). Moreover in the present model for any value of iV 
(N = 2S + 1) we have bound pairs but not complexes of iV bound particles 
as in [16]. 

We conclude the letter with some remarks about the lattice vertex model 
counterpart of the quantum chain considered here. The quantum i?-matrix 
has 1 + 3iV + 2N 2 non-zero Boltzmann weights , which are given by 

R° ° = 1, R° Q 2 = R"o = z sinh A/ sinh( 7 - £l A) 
= Ka = sinh 7/ sinh( 7 - £l A) 
R ~,s = [<W<5/3, 7 + $(i\)u a us6 a , N -p +1 6 SjN - 1+1 \ sinh(7 - A)/ sinh(7 - exA).(21) 
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where a, (3 — 1,2, ...,N. The associated spin Hamiltonian can be found by 
taking the logarithmic derivative of the row-to-row transfer matrix at A = 0. 
It gives the Hamiltonian (3) after a Jordan- Wigner transformation. Since we 
verified that (21) satisfy the Yang-Baxter equations the exact integrability 
of (3) is an immediate consequence. The above vertex model can be treated 
by the diagonal-to-diagonal Bethe ansatz method [29,30], but this is not the 
aim of this letter. 
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